Up and down quark masses from Finite Energy QCD sum rules to five loops 
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The up and down quark masses are determined from an optimized QCD Finite Energy Sum Rule 
(FESR) involving the correlator of axial-vector divergences, to five loop order in Perturbative QCD 
(PQCD), and including leading non-perturbative QCD and higher order quark mass corrections. 
This FESR is designed to reduce considerably the systematic uncertainties arising from the 
(unmeasured) hadronic resonance sector, which in this framework contributes less than 3-4% to 
the quark mass. This is achieved by introducing an integration kernel in the form of a second 
degree polynomial, restricted to vanish at the peak of the two lowest lying resonances. The driving 
hadronic contribution is then the pion pole, with parameters well known from experiment. The 
determination is done in the framework of Contour Improved Perturbation Theory (CIPT), which 
exhibits a very good convergence, leading to a remarkably stable result in the unusually wide 
window so = 1.0 — 4.0 GeV 2 , where so is the radius of the integration contour in the complex energy 
(squared) plane. The results are: m u (Q = 2 GeV) = 2.9 ± 0.2 MeV, m d (Q = 2 GeV) = 5.3 ± 0.4 
MeV, and (m u + m d )/2 = 4.1 ± 0.2 MeV (at a scale Q = 2 GeV). 



The up and down quark masses are very important pa- 
rameters of the QCD sector of the Standard Model. They 
measure the strength of chiral SU(2)®SU(2) and flavour 
SU(2) symmetry breaking, and are believed to be respon- 
sible for the proton-neutron mass difference. For this 
reason, many attempts have been made in the past to 
determine their values in various frameworks, e.g. QCD 
sum rules 0-0, and Lattice QCD 0, while chiral per- 
turbation theory provides information on quark mass ra- 
tios 0. In the framework of QCD sum rules, the ideal 
correlators are those involving the light quark masses as 
overall multiplicative factors, e.g. pseudoscalar correla- 
tors. Unfortunately, the corresponding hadronic reso- 
nance spectral functions are not realistically measurable 
in these channels (e.g. from r decays). Channels where 
they are known from data (vector, axial-vector) involve 
the quark masses as subleading terms, and thus compli- 
cate their determination. The calculation of the scalar 
(pseudoscalar) correlator in Perturbative QCD (PQCD) 
has been improving over the years 0, reaching recently 
the five- loop order [fj . Such level of precison is, sadly, not 
matched in the hadronic sector. In fact, the lack of direct 
experimental data on the resonance spectral functions 
results in a serious systematic uncertainty, which limits 



the accuracy of quark mass determinations in this frame- 
work. While in both the non-strange and the strange 
pseudoscalar channels there are a couple of resonances 
with known masses and widths 0, this is not enough 
to fully reconstruct the spectral functions. In fact, in- 
elasticity and non-resonant background are realistically 
impossible to model. Given this situation, a possible way 
of improving the precision of quark mass determinations 
is to reduce considerably the numerical contribution of 
hadronic resonances. We have recently achieved this in 
the strange pseudoscalar channel 0-0 by using new Fi- 
nite Energy Sum Rules (FESR) involving an integration 
kernel in the form of a second degree polynomial which 
vanishes at the peaks of the two known resonances. This 
resonance quenching has the effect of reducing consider- 
ably (by up to one order of magnitude) the impact of 
this term, so that the pseudoscalar meson pole becomes 
the dominant hadronic contribution. The strange quark- 
mass thus determined is essentially free of systematic 
uncertainties, and subject to future improvement from 
more precise determinations of Aqcd and the gluon con- 
densate, as well as from higher order PQCD calculations 
(six loop and beyond). In addition, the strange quark 
mass at a given fixed scale, e.g. m s (2GeV) turns out 
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to be remarkably stable in a very wide range of values 
of So, the upper limit of integration in the FESR. In 
fact, this range is sq = 1 — 4 GeV 2 , a rather unusually 
broad window. In contrast, without the use of the new 
integration kernel, this window shrinks considerably to 
so = 1 — 2 GeV 2 , and the ghost of the uncontrolable sys- 
tematic uncertainties is resurrected. 
In this paper we use this FESR in the light pseu- 
doscalar channel to find the value of the up and down 
quark masses to five-loop order in PQCD, and includ- 
ing the leading vacuum condensates and higher order 
quark mass corrections. We use the framework of Con- 
tour Improved Perturbation Theory (CIPT) 10|, which 
has proven far more stable than Fixed Order Perturba- 
tion Theory (FOPT) (see 0), at least in the type of 
FESR we are using, with similar conclusions drawn from 
a precision determination of a s (M T ) from r - decay in 
We begin by introducing the correlator of axial- 
vector divergences 

^(q 2 ) =i [d 4 x e*** < \T(d^(x), d v At(0))\ >, (1) 



where d^A^{x) = (ra u + m<j) : d(x) i 75 u(x) : is the 
divergence of the axial- vector current. To simplify the 
notation we shall use in the sequel m u +md = m. Finite 
Energy QCD Sum Rules involving the second derivative 
of this correlator follow from Cauchy's theorem in the 
complex energy-squared, s - plane, Fig. 1 (for details on 
the derivation see [§]) , i.e. 

' ' ds^( S )\ QCD [F(s)-F( So )] 
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where sth = (3iW T ) 2 is the resonance threshold, / ff 
92.4 ± 0.3 MeV, the function F(s) is 
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and the integration kernel ^{s) is 

As(s) = 1 — ao s — ai s 2 

where ao, and a\ are free parameters to be fixed by the 
requirement that A 5 (M 1 2 ) = A 5 (M|) = 0, with Mi j2 the 
masses of the two resonances in the pseudoscalar channel, 
tt(1460) and tt(1830) 0. This gives 

a Q = 0.897 GeV -2 on = -0.1806 GeV -4 
The function -0 5 " (Q 2 ) in PQCD is given by 
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with Q 2 = —q 2 , and Renormalization Group (RG) im- 
provement has been used to dispose of the logarithmic 
terms. The rather long four- and five-loop expressions 
can be found in Since it will turn out that 

so ~ 1 — 4 GeV 2 , we have taken np = 3 in Eq.(6) 
and in the sequel. Using two flavours, instead, changes 
the results by a tiny amount well within the final error 
bars. Contour Improved Perturbation Theory [Io| has 




FIG. 1: Integration contour in the complex s-plane. 



been shown to provide better convergence than FOPT in 
the QCD analysis of the vector and axial-vector correla- 
tors in tau-lepton decay [I3|-[l3|- We found this to be 
also the case in our determination of m s Q, and it will 
be confirmed here for m u ,d ■ Unlike the case of FOPT, 
where a s (so) is frozen in Cauchy's contour integral and 
the (RG) is implemented after integration, in CIPT a s 
is running and the RG is used before integrating. This 
is done through a single-step numerical contour integra- 
tion and using as input the strong coupling obtained by 
solving numerically the RG Equation (RGE) for a s (— s) 
. This technique achieves a partial resummation of the 
higher order logarithmic integrals, and improves the con- 
vergence of the PQCD series. In the case of the pseu- 
doscalar correlator involving the running quark mass as 
an overall multiplicative factor, implementation of CIPT 
requires that not only the running coupling but also the 
running quark mass be integrated around the Cauchy 
contour. The running quark mass can be computed at 
each step by solving numerically the corresponding RGE. 
We define the left hand side of Eq.(2) as 
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(6) and obtain in PQCD (for more details see 0]) 
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FIG. 2: Up quark mass at 2 GeV as a function of so for 
Aqcd = 365 (397) MeV, upper (lower) curve, respectively. 



FIG. 3: Down quark mass at 2 GeV as a function of so for 
Aqcd = 365 (397) MeV, upper (lower) curve, respectively. 



where the RGE for the mass and coupling have been used, 
a s {x) = a s (x)/ir, and where 



F ( x ) = Z^v ) b N s e 

N=l 



(9) 



The constants above are (np = 3): K n = Goi, K\ = 
C'n + 2C\2, K2 — C21 + 2C22, K3 ~ C31 + 2C32, 
K 4 = Go. + 2G 42 , with G01 = 6, G11 = 34, C 12 = -6, 
G21 = -105 C(3) + 9631/24, G 22 = -95, G 23 = 17/2, 
G31 = 4748953/864 - tt 4 /6 - 91519 <(3)/36 + 715 C(5)/2, 
G 32 = -6 [4781/18 - 475 C(3)/8], G 33 = 229, G 34 = 
-221/16, G41 = 33532.26, G 42 = -152 30.6451, G 43 = 
3962.45493, G 44 = -5 34.05 2 83, G 45 = 24.1718750, 
and £(a;) is Riemann's zeta function. Finally, 61 = 
-(s - a So/ 2 - aisg/3), 6 2 = 1/2, fo 3 = -a /6, and 
6 4 = — ai/12. Regarding the value of Aqcd entering 
a s (so); it can be extracted from the strong coupling ob- 
tained from r-decay Q, [HI]. We use the latest high 
precision result [III leading to Aqcd — 365 — 397 MeV. 
The contribution of the gluon condensate to the left hand 
side of Eq.(7) is 
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where the scale /io — 1 GeV 2 appears in connection with 
the removal of logarithmic quark mass singularities (see 
@|). In obtaining this expression use was made of the 
result [H 
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The quark condensate, the dimension six four-quark con- 
densate, and the higher order quark mass corrections 
turn out to be negligible in this application. Turning 
to the hadronic sector, we follow the procedure outlined 
in [14j , where the resonance part of the spectral function 
is written as a linear combination of two Breit-Wigner 
forms normalized at threshold according to chiral per- 
turbation theory, i.e. 
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where BW{s) is a linear combination of two Breit-Wigner 
forms normalized to BW(Q) = 1. The phase space 
factor above is shown in the chiral limit for simplicity. 
The full expression including the threshold at s = 9M 2 
is given in [lj|. Substituting Eqs.(8) and (10) in the 



sum rule Eq.(2), using (%- G )| 



0.06 GeV hi 



Aqcd = 365 — 397 MeV [ll|, and m u /md = 0.553 from 
CPT we find the results shown in Fig. 2 and Fig. 3, 
for m u (2 GeV) and md(2 GeV), respectively. The stabil- 
ity of the results is exceptional, and the duality window 
unusually broad. In order to achieve a reasonable, but 
conservative error, we have analyzed the impact of the 
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various contributions to the quark masses. Some of these 
contributions have opposite effects, e.g. decreasing the 
hadronic resonance term has the same effect as increas- 
ing the QCD contribution, and vice-versa. Allowing for 
a 30% uncertainty (upwards or downwards) in the reso- 
nance spectral function, multiplying (dividing) the gluon 
condensate by a factor of two, and assuming the unknown 
six-loop PQCD term to be equal to the five-loop one, and 
taking into account the current uncertainty in Aqc d , we 
arrive at the following results 

m u {2 GeV) = 2.9 ± 0.2 MeV , 

m d {2 GeV) = 5.3 ± 0.4 MeV 



m ud 



m d 



= 4.1 ± 0.2 MeV 



(13) 



The allowance of a ±30% uncertainty in the hadronic 
resonance contribution impacts the quark masses at the 
level of only 1 %. Our choice of integration kernel, Eq. (4) , 
is by no means unique. We have tried several other func- 
tional forms with the same constraint, in particular the 
so called pinched kernels 15]. These kernels are problem- 
atic in this application, as the relative contribution of the 
poorly known spectral function is much larger, and the 
stability region is considerably reduced. The same effect 
was observed in the case of the strange quark mass [9]. 
We have also checked that using polynomials of higher 
order, again restricted to vanish at the resonance peaks, 
has a negative impact on the QCD sector of the FESR 



due to the emergence of unknown higher dimensional vac- 
uum condensates, and the divergence of polynomial coef- 
ficients. From the point of view of errors our polynomial 
parametrization, Eq. (4), represents the optimal choice. 
Disregarding the stability region as a criterion, results 
for the quark masses with other choices for the kernel 
are consistent within errors with the results above. In 
particular, for a given hadronic resonance spectral func- 
tion, e.g. Eq.(12), if we remove the kernel altogether 
(A 5 (s) = 1) the stability region shrinks to the narrow 
window so = 1 — 2 GeV 2 , and the systematic uncertain- 
ties, from the (unmeasured) hadronic resonance spectral 
function, make a comeback. 

In comparing these results with others from QCD sum 
rules, one should exercise extreme care, as older determi- 
nations used much lower values of Aqc d , and less terms 
in PQCD (with the exception of the first five-loop re- 
sult of Q). However, the most important difference with 
previous determinations lies in our choice of integration 
kernel, which reduces considerably systematic uncertain- 
ties from the resonance sector. The results above com- 
pare favourably with some recent lattice QCD determina- 
tions, e.g. for the average: m u d = 3.85 ± 0.42 MeV [lq . 
jtw = 3.72 ± 0.41 MeV [13], and m ud = 3.3 ± 0.2 MeV 
[18j, where all errors have been combined in quadra- 
ture. However, one other lattice QCD result of m u d = 
2.527 ± 0.047 MeV is much smaller than these. 
This work was supported in part by NRF (South Africa) 
and DFG (Germany) 



[1] For a recent review see e.g. P. Colangelo, A. Khod- 
jamirian, in: "At the Frontier of Particle Physics/ Hand- 
book of QCD'", M. Shifman, ed. (World Scientific, Sin- 
gapore 2001), Vol. 3, 1495-1576. Some recent papers in- 
clude: S. Narison, Phys. Rev. D 74, 034013 (2006); E. [6] 
Gamiz, M. Jamin, A. Pich, J. Prades, F. Schwab, Phys. 
Rev. Lett. 94, 011803 (2005) ; M. Jamin, J.A. Oiler, A. 
Pich, Phys. Rev. D 74, 074009 (2006), Eur. Phys. J. C [7] 
24, 237 (2002) ; K. Maltman, J. Kambor, Phys. Rev. D 
65, 074013 (2002). [8] 

[2] M. Jamin, M. Miinz, Z. Phys., Particles & Fields C 66, 

633(1995); K.G. Chetyrkin, OA. Dominguez, D. Pirjol, [9] 

K. Schilcher, Phys. Rev. D 51, 5090 (1995); K. Chetyrkin, 

D. Pirjol, K. Schilcher, Phys. Lett. B 404, 337 (1997). [10] 

[3] HPQCD Collaboration, C. Aubin et al., Phys. Rev. D 70, 
031504 (2004); CP-PACS and JLQCD Collaborations, T. 
Ishikawa et al., arXiv: |hep-lat/050"9"T42~| PoS LAT 2005, 
057 (2006); QCDSF Collaboration, M. Gockeler at al., [11] 
arXiv: |hep-lat/0509~i"59| PoS LAT 2005, 078 (2006); D. 
Becirevic et al., Nucl. Phys. B 734, 138 (2006); Q. Mason, [12] 
H.D. Trottier, R. Horgan, C.T.H. Davies, G.P. Lepage, 
Phys. Rev. D 73, 114501 (2006). [13] 

[4] H. Leutwyler, Phys. Lett. B 378, 313 (1996). 

[5] K.G Chetyrkin, A.L. Kataev, F.V. Tkachov, Phys. Lett. [14] 
B 85, 277 (1979); M. Dine.J. Sapirstein, Phys. Rev. Lett. 
43, 668 (1979) ; W. Celmaster, R. Gonsalves ibid. 44, 560 



(1980); S.G Gorishny, A.L. Kataev, S.A. Larin, Phys. 
Lett. B 259, 144 (1991) ; LR. Surguladze, M. Samuel, 
Phys. Rev. Lett. 66, 560 (1991); T. van Ritbergen, J.A.M. 
Vermarseren, S.A. Larin, Phys. Lett. B 400, 379 (1997). 
PA. Baikov, K.G. Chetyrkin, J.H. Kuhn, Phys. Rev. 
Lett. 96, 012003 (2006); K.G. Chetyrkin, A. Khod- 
jamirian, Eur. Phys. J. C 46, 721 (2006). 
Review of Particle Physics, Particle Data Group, J. Phys. 
G: Nucl. Part. Phys. 33, 1 (2006). 

OA. Dominguez, N. Nasrallah, K. Schilcher, J. High En- 
ergy Phys. 0802, 072 (2008) . 

OA. Dominguez, N. Nasrallah, R. Rontsch, K. Schilcher, 

J. High Energy Phys. 0805, 020 (2008). 

A. A. Pivovarov, Z. Phys., Particles & Fields, C 53, 461 

(1992); F. Le Diberder and A. Pich, Phys. Lett. B 286, 

147 (1992); M. Jamin, J. High Energy Phys. 0509, 058 

(2005). 

M. Davier, S. Descotes-Genon, A. Hocker, B. Malaescu, 
Z. Zhang, arXiv: 0803.0979. 

OA. Dominguez, K. Schilcher, J. High Energy Phys. 
0701, 093 (2007). 

M. Davier, A. Hocker, Z. Zhang, Rev. Mod. Phys. 78, 
1043 (2006). 

OA. Dominguez, Z. Phys., Particles & Fields, C26, 269 
(1984), OA. Dominguez, E. de Rafael, Ann. Phys. (New 
York) 174, 372 (1987). A misprint in the latter has been 



■5 



corrected in: J. Bijnens, J. Prades, E. de Rafael, Phys. 
Lett. B 348, 226 (1995). 
[15] K. Maltman, Phys. Lett. B 440 (1998) 367; C.A. 
Dominguez, K. Schilcher, Physics Letters B 448, 93 
(1999); ibid. B 581, 193 (2004). 



[16] B. Blossier et al., J. High Energy Phys. 0804, 020 (2008). 

[17] C. Allton et al., arXiv: 0804.0473. 

[18] C. Bernard et al. PoSLAT 2007, 090 (2007). 

[19] S. Aoki et al., arXiv: 0807.1661 



